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I. INTRODUCTION 



The recent discoveries of the E&, and baryons at the Tevatron have incited to a 
new analysis of heavy baryons both within the combined 1/N C and 1/rriQ expansion jl| and 
the quark model, see for example Refs. . The combined 1/N C and 1/rriQ expansion 

is a model independent method. It is thus important to search for a link between this 
method and the quark model. In previous studies [5,^] we have investigated the possibility 



to establish a connection between the two approaches and we have : 



bund that a remarkab 



compatibility exists between them when dealing with nonstrange 5j or strange baryons 6]. 

Presently we extend the ideas of our previous studies 5, (| to the case of heavy baryons 
made of one heavy quark (c or b) and two light ones (u, d, or s). This is the first step of a 
larger project and we view it as an exploratory work where we search for the compatibility 
between the spin-independent part of a quark model Hamiltonian and the corresponding 
terms in the combined 1/N C and 1/m.g expansion mass formula for the ground state. The 
spin-dependent part as well as the excited states will be analyzed subsequently. 

As previously, the comparison of the quark model results with those of the 1/N C expansion, 
presently combined with an 1/ vtlq expansion, will be based on the introduction of a quantum 
number N, which is the same as in the harmonic oscillator potential and which is treated as 
a band number in baryon phenomenology. The introduction of N in the eigenvalues of the 
Hamiltonian was quite simple for identical quarks, the procedure becomes more involved for 
baryons containing heavy quarks, as we shall see. 

The paper is organized as follows. After a summary of the charm and bottom baryon 
flavor states given in Sec. [Ill the mass formula used by combining the 1/N C and 1/tuq 
expansions for such baryons is presented in Sec. IIHI Section [IV] gives a corresponding 
mass formula obtained from a Hamiltonian quark model where the confinement is of Y- 
junction type and where one gluon exchange and quark self-energy contributions are added 
perturbatively. In that section the excitation quantum number N is introduced and its 
meaning is discussed. A comparison between results obtained on one hand in the combined 
1/rriQ and 1/N C expansion and on the other hand in the quark model is then made in Sec. [V] 
Conclusions are finally drawn in Sec. IVI1 

In the following, the symbol q will denote a light quark (u, d, s) and the symbol Q will 
denote a heavy quark (c, b). Moreover, the symbol n will be used for u and d quarks since 
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TABLE I: Flavor states of the single charm sextet baryons: J p = | (no star) for baryons in the 
mixed representation and J p = | + (with star) for baryons in the symmetric representation. They 
all have J qq = 1. Members of the same doublet become degenerate at uiq —* oo. 

Baryon doublet Flavor state 

£++, £*++ uuc 

S+, £*+ -^{ud + du)c 

E°, £*° ddc 



(us + su) c 
(ds + sd) c 



0°, ft* ssc 



both particles are assumed to have the same mass, as in our previous works. 



II. FLAVOR STATES 
A. Charm baryons 

Here we introduce the classification of ground state heavy baryons based on SU(4). In 
the following the total spin of a baryon is denoted by J, the spin of the light subsystem 
by Jqq, and that of the heavy quark by Jq. In SU(4) the baryon multiplets arise from the 
direct product decomposition 4 © 4 © 4 = 20 © 20 © 20 © 4, see e. g. Ref. [i| . All baryons 
in the symmetric multiplet 20 have J p = | + . The lightest SU(3) submultiplet is the well 
known Gell-Mann-Ne'eman decuplet. The single charm baryons form a sextet where the 
Fermi statistics requires J qq = 1. The six baryons ^* c ++ , "*' + ; ^* ° an d have 

the flavor structure given in Table [B The remaining members of the symmetric multiplet 
are the three double charm baryons and the triple charm baryon ^++ + . 

The experimental masses of single charm baryons with J p = | + are 

£* = 2518.0 ± 0.8 MeV, 
S* = 2646.4 ± 0.9 MeV, 

tt* c = 2768.3 ± 3.0 MeV. (1) 
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TABLE II: Flavor states of the single charm antitriplet baryons with J 1 



| in the mixed 



symmetric representation. They all have J qq = 



Baryon 



Flavor state 



A c + 



\ [ud — du) c 
\ (us — su) c 
i (ds — sd) c 



which represent mass averages when the hadron appears with different charges. Note that 
here and below none of the quantum numbers assigned to the charm baryons have been 
measured experimentally, but are based on quark model expectations. 

The mixed symmetric representation 20 has J p = | + . The lowest submultiplet is the 
SU(3) Gell-Mann-Ne'eman octet. The single charm baryons £+, E^ + , S^ , and f2° 
form a sextet with flavor states indicated in Table [land J qq = 1. A+, H+, and S° form an 
antitriplet with flavor states shown in Table [III having J qq = 0. 

The experimental masses of single charm baryon with J p = | + are 

A c = 2286.46 ±0.14 MeV, 

S c = 2453.56 ±0.16 MeV, 

E c = 2469.5 ± 0.3 MeV, 

E' c = 2576.9 ±2.1 MeV, 

Q c = 2697.5 ± 2.6 MeV, (2) 

where, again, mass averages are made when the hadron appears with different charges. In 
the observed S c and H' c baryons it is expected that the light quarks are mostly in a state 
with J qq = and J qq = 1 respectively. 

The mixed symmetric multiplet also contains three double charm baryons H+, and 
f2+ from which only 5+ has been observed by SELEX with a mass of 3518. 9±0. 9 MeV Q], 
but needs confirmation. 
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TABLE III: Flavor states of the single bottom sextet baryons: J p = ^ + (no star) for baryons in 
the mixed representation and J p = | + (with star) for baryons in the symmetric representation. 
They all have J qq = 1. Members of the same doublet become degenerate at mg — > oo. 

Baryon doublet Flavor state 

E^ , uub 
Eg, E*° ^( ud + du ) h 

^ (us + su) b 

3r> S 6*'~ -±(ds + sd)b 

O^, ssb 

TABLE IV: Flavor states of the single bottom antitriplet baryons with J p = ^ + in the mixed 
symmetric representation. They all have J qq = 0. 

Baryon Flavor state 



, E* b ~ ddb 
=•'0 ^*'o 1 

' ^b 



A^ <J \ (ud - du) b 

| (us — su) b 
fx 

-6 V 2 



-° '± (ds - sd) 6 



B. Bottom baryons 

Despite the large symmetry breaking, for the sake of the classification one can also assume 
an SU(4) classification of bottom baryons. Similarly, for single bottom baryons there is a 
sextet shown in Tabl eJTTl and an antitriplet shown in Table IIVI The mass of A& has been 
previously measured [8| 

A b = 5620.2 ± 1.6 MeV. (3) 
Recent measurements have been made for H& 0, Q, E&, and U|. The measured 



5 



masses are 



5774 ± 11 ± 15 MeV 



5811.5 ± 1.7 MeV 



5832.7 ± 1.8 MeV [11]. 
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% 5792.9 ± 2.5 ± 1.7 MeV [loj] , 



(4) 



The remaining undiscovered single bottom baryons are S' b , H£ , fit, and 



III. GROUND STATE HEAVY BARYONS IN THE l/m Q AND 1/N C EXPANSION 

The approximate spin-flavor symmetry for large iV c baryons containing light q = u,d, or s 
quarks and heavy Q = c or b quarks is SU(6)x SU(2) C x SU(2)&, i.e. there is a separate spin 
symmetry for each heavy flavor. Over a decade ago the 1/N C expansion has been generalized 



3. 



to include an expansion in 1/rriQ and light quark flavor symmetry breaking 

Let us first consider that SU(3)-flavor symmetry is exact. In this case the mass operator 
is a flavor singlet. In the combined I/tuq and 1/N C expansion to order I/tUq the ground 
state mass operator takes the following form 

M (1) = m Q N Q l + A qq + A Q + X qqQ , (5) 

where Nq is the number of heavy quarks. The leading order term is rag at all orders in the 
1/N C expansion. Next we have 



A qq = Co N c t + ^-J 2 qq , (6) 



where J qq is the total spin of the light quark pair. This operator contains the dynamical 
contribution of the light quarks and is independent of itlq. Then, \q gives the I/tuq 
corrections due to Nq heavy quarks 

In the following we shall deal with Nq = 1 only. Lastly, X qq Q contains the heavy-quark 
spin-symmetry violating (chromomagnetic) operator which is of order I/tuq as well 

// 

\qQ = 2 AT — J qq ' Jqi ( 8 ) 
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TABLE V: Mass combinations resulting from heavy quark and large N c limit and their experimental 
values 



Mass combination 


Exp (MeV) 


Exp (MeV) 




Q = c 


Q = b 


Aq 


2286.46±0.14 


5620. 2±1. 6 


i(E Q + 2E^)-A Q 


210.0±0.5 


205.4±2.1 




64.4±0.8 


21.2±2.5 




183.0±0.3 


172.7±3.4 



Jq being the spin of the heavy quark. This is the term responsible for the splitting between 
states which form degenerate doublets in the heavy quark limit, see Tables HI and IIHI 

The unknown coefficients Co, C2, c' Q , c' 2 , and c 2 ' are functions of 1/N C and of a QCD scale A. 
Each coefficient has an expansion in 1/N C where the leading term is of order unity and does 
not depend on 1/ttiq. Without loss of generality one can set Cq = A. The other coefficients 
contain a dimensional power of A and a dimensionless function of 1/N C beginning at order 



unity and have to be fitted to the available experimental data. In agreement with Ref. 
we can take 



it 



c = A, c 2 ~ A, 

C ~ C 2 ~ C 2 ~ A ' ( 9 ) 



At the dominant order, the value of A can be extracted from the mass combinations 

A Q = m Q + N C A, (10a) 
i(£ Q + 2^)-A Q = 2A (10b) 

= (1 ° C) 
resulting from the mass definition (jSj). The equations (jlOaj) and (llObj) express the fact that 



Xq is negligible with respect to the other terms in (jSJ). Here and below the particle label 
represents its mass. 

A slightly more complicated mass combination, involving light baryons as well as heavy 
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ones, directly leads to rriQ, that is 



in 



-(2iV + 3E + A + 2S) 



— (4A + 3£* + 22* + ft) 



m Q- 



This mass combination gives 



while the value 



m c = 1315.1 ± 0.2 MeV, m h = 4641.9 ± 2.1 MeV, 



A w 324 MeV 



(12a) 



(12b) 



ensures that the mass combinations ffTUl) are optimally compatible with the experimental 
values for Q = c and Q = b indicated in Table |Vj Note also that the heavy quark flavor 
symmetry predicts that the observed (A& — A c ) = 3333.7 ±1.6 MeV splitting [l| can give a 
measure of the quark mass difference m& — m c up to corrections of the order A 2 (l/2m c — 



l/2m&) ~ 23 MeV 121 ] . The values given by Eqs. ( 1121) satisfy this constraint. 

The operator analysis including SU(3)-flavor breaking leads to an expansion in the SU(3) 
violating parameter e which contains the singlet M' 1 ', an octet and a 27-plet M^ 27 \ 

The last term brings contributions proportional to e 2 and we neg lect it. For we retain 
its dominant contribution T 8 to order N®. Then the mass formula becomes 



M = M (1) + eT 8 . 



(13) 



The flavor breaking parameter e is governed by the mass difference m s — m (where m is the 
average of the m u and masses) and therefore is e ~ 0.2-0.3. It is measured in units of 
the chiral symmetry breaking scale parameter A x ~ 1 GeV. A measure of the SU(3)-flavor 



breaking factor can be given by 12J 



Ho- A 



Q 



(14) 



The value (eA x ) = 206 MeV leads to Sq — Aq = 178 MeV, which is the average value of the 
corresponding experimental data listed in Table |V] 

IV. QUARK MODEL FOR HEAVY BARYONS 
A. Hamiltonian 



The potential model used to describe heavy baryons is the same as that which has been 



proposed in Ref. 



6] for light baryons. Let us recall its main features. 
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In quark models, a baryon is a bound state of three valence quarks which can be described 
at the dominant order by the spinless Salpeter Hamiltonian 

3 



i=i 

where m 8 is the current (bare) mass of the quark i and Vy is the confining interaction 
potential. Both the flux tube model [3| and lattice QCD Q suggest that the flux 
tubes form a Y-junction: A flux tube starts from each quark and the three tubes meet at 
the Torricelli point of the triangle formed by the three quarks. This point, located in xt, 
minimizes the sum of the flux tube lengths, leading to the following confining potential 

3 

Vy = a^^\xi — xt\ ■ (16) 

i=i 

The position of the quark i is denoted by Xi, and a is the energy density of the flux tubes. 
Such a Hamiltonian can also be obtained in the framework of the field correlator method 

As Xt is a complicated three-body function, it is interesting to approximate the confining 
potential by a more tractable form. In the following, we will use 

3 



-i=i 

3 

V R = ka^2\xi-R , (18) 



where R is the position of the center of mass and k is a corrective factor [17|. The eigenvalues 
corresponding to potentials Vy and Vr differ from each other only by about 5% in most cases. 
The accuracy of the formula (ITS]) is thus rather satisfactory, and has already led to relevant 
results in Ref. jf|. For light (symmetrical) qqq baryons, a good value for the corrective factor 
is ko = 0.952. For very asymmetrical qqQ baryons, a good choice is k\ = 0.930 [13]. This 
last value corresponds actually to the case m q /rriQ — ► 0. 

Besides the confining potential ({TBI) , other contributions are necessary to reproduce the 
baryon masses. We shall add them as perturbations to the dominant Hamiltonian (fT7|) . The 
most widespread correction is a Coulomb interaction term of the form 



arising from one gluon exchange processes, where as^j is the strong coupling constant be- 
tween the quarks i and j. Actually, one should deal with a running form as(r), but it would 
considerably increase the difficulty of the computations. Typically, we need two values: 
«o = a s,qq for a 11 pair and a\ = cts,qQ for a qQ pair, in the spirit of what has been done 
in a previous study describing mesons in the relativistic flux tube model [3]. There it was 
found that ai/a ~ 0.7 describes rather well the experimental data of qq and Qq mesons. 

Another perturbative contribution to the mass is the quark self-energy. This is due to 
the color magnetic moment of a quark propagating through the vacuum background field. It 



adds a negative contribution to the hadron masses 191 ] . The quark self-energy contribution 
for a baryon is given by 

The factors / and 6 have been computed in quenched and unquenched lattice QCD studies 



m - Although a 15 ,lot k,lown - rf eat accuracy - a seems we " established that 



3 < / < 4 and (1.0 < S < 1.3) GeV [2fJ, ]2l|. The function 77(e) is analytically known; we 
refer the reader to Ref. [3] for the explicit formula. For typical values of the light quark 
masses, we have < m q /5 < 0.3, while for heavy quarks, we have 1.0 < rriQ/5 < 6.0. The 
function 77(e) is such that 



rj(e) w 1+ (4 + 3 In e 2 for e < 1, 



2 
72 



— for 00. (21) 



e 

For the relevant values of e = rrij/5 a better accuracy is obtained with the following simple 
forms 

77(e) w 1 - /5e 2 with /3 = 2.85 for < e < 0.3, 



7 

c 2 



~ — with 7 = 0.79 for 1.0 < e < 6.0. (22) 



e z 

3 



Let us note that the corrections depending on the parameter 7 appear at order 1/ttiq in the 
mass formula, so they are not considered in this work. Finally, fii is the dynamical mass of 
the quark i, defined as pjj] 



^ = \\jp! + ™ij- (23) 

This dynamical mass is state-dependent: It represents the kinetic energy of the quark i 
averaged with the wave function of the unperturbed spinless Salpeter Hamiltonian ffTTl) . 
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B. General formulas 



We are mainly interested in analytical expressions, so that a comparison with the large 



N c mass formula will be straightforward. To this aim, the auxiliary field technique wil 



used in order to transform the Hamiltonian (1171) into an analytically solvable one 22 
With A = k a, we obtain 



be 



23] 



3 

E 



p f + m 2 



2u 3 



3 



\ 2 (xj - R) 2 Vj 
2^ + 2~ 



(24) 



The auxiliary fields, denoted as /ij and Vj, are operators, and H(fii, Vj) is equivalent to 
up to their elimination thanks to the constraints 








p? + m h 



\\xj — R\. 



(25) 



(fJiifl) is the dynamical quark mass introduced in Eq. ( l23i) . and (^,0) is the energy of the flux 
tube linking the quark i to the center of mass. 

Although the auxiliary fields are operators, the calculations are considerably simplified 
if one considers them as real numbers. They are finally eliminated by a minimization of 
the masses 22], and the extremal values of /ij and Vj are logically close to (/Xj,o) and (^,0) 
respectively. This technique can give approximate results very close to the exact ones (see 
Ref. 24j for a comparative study of baryons with the auxiliary fields introduced only in the 
kinetic part of the Hamiltonian) . 



In Ref. 



251 ] , it has been shown that the eigenvalues of a Hamiltonian of the form (|24p can 



be analytically found by making an appropriate change of variables, the quark coordinates 
Xi = {^1,^2,^3} being replaced by new coordinates x k = {R,£,ff}. The center of mass is 
defined as 

R= , (26) 



/it 

— # — # 

with [i t = 1^2+1^3- {£, V } are two relative coordinates: £ cx X1—X2 and 77 cx 



-x 3 . 



As we only consider baryons built from two different quarks, the general formulas obtained 



in Ref. 25j can be simplified. In the case of two quarks with mass m and another with mass 
777,3, the mass spectrum of the Hamiltonian (124]) is given by (/ii = /i 2 = /J, V\ = ^2 = u ) 



M(/i, /i 3 , v 3 ) = uz(Nz + 3/2) + u v (N v + 3/2) + \l + 1/ + 



A*3 + ^3 



in 



/i 2^3 



(27) 
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= — , , X J^ + ^- (28) 



where 



The integers N^/ v are given by 2n^/ v + £^/ v , where n^/ v and £^/ v are the radial and orbital 
quantum numbers relative to the variable £/ff respectively. One can also easily check that 
25 1 

?)^, (29) 
with 



(30) 

2/i + /i 3 

These last identities provide relevant informations about the structure of the baryons, since 
<*'> = «*-*>'>=V=g>>- (3D 

Moreover, by symmetry, we can assume the following equality 

<(* a - f 3 ) 2 > = {{x 2 - f 3 ) 2 > « ^ + (Y 2 ) , (33) 

which will be useful in the computation of the one gluon exchange contribution. 

The case of qqq baryons, studied in our previous papers , is obtained by taking 

m = m n = and m 3 = m s , and by setting A = k a. If the three quarks are identical, 
then m 3 = m, /i 3 = /z, 1/3 = ZA For qqQ baryons, we explicitly write m 3 = ttiq, /x 3 = /iq, 
^ 3 = ^Qj and we set A = k\a, m = or m s for n or s quarks respectively. Let us note that 
different values of ko have been previously used: ko = (1/2 + \/3/4) in Ref. [5| and k = 1 
in Ref. 6|. In this work, we choose phenomenological values computed in Ref. 17| in order 
to obtain the best possible simulation of the Y-junction for both qqq and qqQ baryons with 
the potential fTTSj) . 



C. Mass formula for heavy baryons 

In this section, we focus our attention on ssQ baryons. The mass formula for nnQ 
baryons is obtained simply by setting m s = 0, and the case of nsQ baryons will be discussed 
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in the next section. The four auxiliary fields appearing in the mass formula (12 7p have to be 
eliminated by solving simultaneously the four constraints: 



d^M(fi,fiQ,u,u Q ) = 0, d^ Q M(fi,fi Q ,u,u Q ) = 0, 

d v M(ii,iiQ,v,v Q ) = 0, d UQ M(ji,fiQ,u,v Q ) = 0. (34) 

This cannot be done exactly in an analytical way, but solutions can be obtained by assuming 
that 1 / tuq and m s are small quantities. After some algebra, a solution was found by working 
at order I/thq and m 2 s (all contributions proportional to m s are vanishing). By denoting 

N = N^ + N v , 



'kMN + 3) 

Mi 

G(N, N„) = y/2N v + 3 (v/2(iV + 3) - y/2N v + 3) , (35) 

we have obtained 

3m 2 „ kia Ar Ar . 

„ = «i _ *i!L (2iv + 3) 

4/ii 4m Q 
MQ = ^Q + 7T 1 ^G'(iV,iV ?? ), 



vq = ^l(2N n + 3)(N + 3) (3(J) 

Logically, /ig w since this auxiliary field is dominated by the effective mass of the heavy 
quark. The length of the flux tube joining the heavy quark to the center of mass is smaller 
than the other ones, so linim^oo vq = as expected. 

The mass formula (1271) . in which the auxiliary fields are replaced by the expressions (1361) . 
reads at orders 1 / ttiq and m 2 s as 

2 1 

M = m Q + 4/i! + ^ + -^G(iV, N v ). (37) 
/ii 2mQ 

It is interesting to look at the magnitude of the various terms in this formula. Let us choose 
typical values for the parameters: ki — 1, a — 0.2 GeV 2 , m s = 0.3 GeV, m c = 1.5 GeV, 
nib — 5.0 GeV. For the ground state (N = 0), /ix = 0.548 GeV. The contribution of the 
kinetic energy and of the confinement in M, given by Afii = 2.191 GeV, is of the order of tuq. 
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The contribution of the strange quark is given by ^ = 0.164 GeV, while the term ^f-G(0, 0) 
is 0.083 GeV and 0.025 GeV respectively for the charm and bottom masses. These values 
justify a posteriori the use of the power expansion in m s and in I/tuq. Formulas (1361) and 
( 1371) giving the optimal values of the auxiliary fields and the corresponding minimal mass are 
approximate solutions of Eq. (T27J). In Table IVa these values are compared with the exact 
solutions obtained numerically. In all cases, the error on the mass is quite small, even if the 
error on some auxiliary fields is larger. The auxiliary fields \x and \xq are used to compute 
perturbatively the self-energy. Fortunately, the error on these fields are small. As expected, 
the accuracy is improved for large values of rriQ, while m s has only a little influence. 



TABLE VI: Relative error (%) on auxiliary fields (|36|) and mass (|37|) for typical values of the 
physical parameters {k\ = 1, a = 0.2 GeV 2 ). Quark masses are given in GeV. 



m s 1 m Q 




/ 1.5 






0.3 / 1.5 






/ 5.0 






0.3 / 5.0 




(N,N V ) 


(0,0) 


(4,0) 


(4,4) 


(0,0) 


(4,0) 


(4,4) 


(0,0) 


(4,0) 


(4,4) 


(0,0) 


(4,0) 


(4,4) 




0.007 


6.6 


5.2 


2.8 


5.6 


5.4 


0.2 


0.7 


1.4 


2.9 


0.02 


2.0 


V 


8.9 


7.1 


29.5 


10.2 


7.2 


30.1 


1.1 


1.0 


3.6 


2.3 


1.1 


3.5 




0.2 


4.8 


4.6 


0.8 


4.5 


5.3 


0.04 


0.2 


0.4 


0.1 


0.1 


0.5 


V Q 


44.5 


74.0 


67.4 


33.7 


68.8 


61.9 


12.8 


21.8 


18.7 


1.8 


16.4 


13.3 


M 


0.2 


1.4 


1.7 


0.5 


1.7 


1.7 


0.05 


0.1 


0.3 


0.6 


0.3 


0.2 



The contribution of the one gluon exchange term can be computed with the help of 
relations ( !31~1) and ( !32|) . One obtains 



AM oge « 




2«! 



+ Y 



4/i| 8m Q 



'kia 



%a 2N v + 3 



4/x 2 2m Q ^2(V + 3) 



'2(2N V + 3) 
V + 3 



- 1 



(38) 



For values of the parameters defined above together with cto = ct\ = 0.4, the contribution of 
the dominant term in AM oge is —0.264 GeV for the ground state. The m 2 s / /if term brings 
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—0.020 GeV while the l/m c term brings 0.034 GeV. Again, the use of the power expansion 
in m s and in l/mq seems relevant. 

The relations ( 1361) defining \i and \iq allow to write down the contribution of the quark 
self-energy ( 1201) . Using the approximation (122"]) one obtains 



AM qse = — 



1 ( A + ll ml + -^G(N,N v ) 
1 4/if 5 2 y 4/xim Q VJ 



(39) 



7T/il 

We recall that the correction proportional to (3m 2 comes from a convenient parameterization 
of the r](e) function, while the term proportional to m 2 / fj, 2 is due to the expansion of the 
auxiliary field \i. For the values of the parameters defined above and the typical values 
/ = 3.5 and 5 = 1 GeV, the contribution of the dominant term in AM qse is —0.302 GeV for 
the ground state. The m 2 / term brings 0.092 GeV while the l/m c term brings —0.031 GeV. 
The use of the power expansion in m s and in I/tuq seems here more questionable, mostly for 
the contribution of the strange quark. This is due to the par ticular nature of the self-energy 
interaction which can be defined only as a perturbation 19]. 

If we now look at the dominant terms in M — ttlq, AM oge and AM qse , we find respectively 
2.191 GeV, —0.264 GeV and —0.302 GeV for the ground state with parameters defined above. 
These numbers show that it is a posteriori justified to treat the Coulomb interaction and 
the self-energy interaction as perturbations. 



D. Mass formulas for general qqq and qqQ baryons 



In this section we gather mass formulas obtained for both light and heavy baryons. The 
qqq mass formula is given in Ref. f| and is reminded here for completeness 
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qqq 
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(40) 



All parameters were already presented above, except the number n s of s-quarks in the 
baryons. The mass formula M qqq depends only on N = + N v since the contribution of 
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terms proportional to Ng — N v , vanishing for n s = and 3, was found to be very weak in 



general 



6] 



In the previous section, only the case of a heavy baryon containing two identical light 
quarks was treated (n s = or n s = 2). It has been shown that every s quark brings the same 
contribution AMq s to the mass of a light baryon [see Eq. (jffij)]. So, we can reasonably assume 
that the same situation occurs for qqQ baryons. To take into account the contribution of 
n s quarks to the mass of these baryons, it is enough to replace the term m? s by n s m 2 s /2 in 
Eqs. (1371) . (1351) and (I3U]) . Let us note that it is not necessarily true for the auxiliary fields \x 
and v [6j. In the following formulas, we keep explicitly the dependence on both Ng and N v : 
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(41) 



E. What is the good quantum number? 



At the lowest order, the mass formula (1371) . with the rescaling a <-> a (see next section), 
leads to 



(M - m Q f 



47T(T fci 

3 fc n 



(iV + 3). 



(42) 



The model thus predicts Regge trajectories for heavy baryons, with a slope of 4naki/ (3k ) w 
1.37TO" instead of 27ra for light baryons. At this dominant order, the mass formula depends 
only on N. However, when corrections are added, the mass formula is no more symmetric in 
N v and N^. Is it still possible to find a single quantum number? There are three possibilities: 

• As in Ref. [6_|, we could assume that N% w N v . But, the presence of a heavy quark 
makes the system rather asymmetric in the £ and fj variables. So this solution seems 
unnatural. 
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Another possibility is to impose N v = and = N. With no excitation in the fj vari- 
able, the two light quarks are moving around a static heavy quark in the configuration 



q — Q — q, as proposed in Ref. [26]. 



• The opposite possibility can also be assumed: N v = N and N% = 0. With no excitation 
in the £ variable, the two light quarks behave as a diquark orbiting around the heavy 
quark by forming a Q — (qq) system, as considered in Ref. [27|. 

At order I/tuq, the dominant term 0371) depends on the function G(N,N V ). The baryon 
mass is lowered when G(N, N v ) is minimal, that is to say for N v = N. In this case 



F(N) = G(N, N) = V3 + 2iV ^2(3 + N) - + , (43) 

with F(0) = 3( v / 2 - 1) < F(N) < 3/2, this upper bound being the limit of F(N) for 
going to infinity. The analysis of the dominant part of the Coulomb term (|38il shows that 
the baryon mass is also lowered when N v = N. So it is natural to assume that the favored 
configuration, minimizing the baryon energy, is N v — N and N% = 0, as in Ref. [27]. In this 
case a light diquark-heavy quark structure for the baryon is favored. 

It is also possible to reach the same conclusion by looking at the mean values of the 
variables X and Y. At the dominant order, we have 

(^-^- 

Because of the particular shape of the potential (a Cornell type), the more the system is 
small, the more its mass will be small. Indeed, the energy of the flux tubes increases with 
the size of the baryon, while the attractive Coulomb-like forces are larger for small quark 
separations. Equations (jUj) shows that an excitation of type N v will keep the baryon smaller 
than the corresponding excitation in N%. Thus, the most favored possibility, at least for the 
small excitation numbers, is also N v = N and N% = 0. 

As for light baryons, heavy baryons can be labeled by a single harmonic oscillator ex- 
citation number and the emergence of this quantum number can be understood within a 
relativistic quark model framework. However, we only discuss the ground state in the fol- 
lowing, that is N% — and N v — N — 0. Excited states will be studied in subsequent 
papers. 
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TABLE VII: Parameters for qqq baryons. 



Fixed parameters Fitted parameters 

m n = a = 0.4 m s = 0.240 GeV 

k = 0.952 5 = 1.0 GeV a = 0.165 GeV 2 

a = 7ro-/(6fc ) /? = 2.85 / = 3.60 



F. Determination of the parameters 

The parameters needed for qqq baryons have been obtained in our previous papers j^, [f| 
but, since we use a new value for ko, we prefer to determine a set of new values for the 
parameters which are gathered in Table IVHI The new values are very close to the previous 
ones and do not alter the good results obtained in Refs. 

Bfl 

. The auxiliary field method 



systematically overestimates the absolute scale of the mass spectrum [22] . In order to obtain 

a good accuracy for the baryon masses, it is necessary to perform the rescaling a = 7r<r/(6fco) 

Pi 

throughout the mass formulas, where a is the physical string tension for a meson [5!]. As 
u and d current quark masses are expected to be very small, we also take a vanishing 
current mass for the quark n. The parameter a and / are fitted on the nnn baryon Regge 



or a a 
5|. It is 



trajectory. As it is not possible to determine independently a and /, we choose 
value in agreement with other potential models. More details can be found in Ref. 
worth noting that the value 3.6 for / is in the range [3-4] and that the string tension value 
of 0.165 GeV 2 is in good agreement with the value predicted by the flux tube model 28]. 
The s-quark mass is fitted to the strange baryon masses in the band N = [6]. The value 
found for m s is larger than the PDG value of 104^ 2 4 MeV [8|. However, a strange quark 



mass in the range 0.2-0.3 GeV is quite common in potential models 29_|, |30|, |31|. 



The parameters linked to heavy quarks are m c , rrib, ki, and a±. We fix a.\ = 0.7ao from 



the quark model study of Ref. [LSj. The value k\ = 0.930 has been computed in Ref. [IT ] . 
Because of the rescaling a = ira/(6ko), only the ratio ki/ko ~ 0.98 is relevant. Let us note 
that to fix this ratio to 1 does not change noticeably the other parameters. The heavy quark 
masses can be fitted to the experimental data as follows. The quark model mass formula (141 p 
is spin- independent; it should thus be suitable to reproduce the masses of heavy baryons for 
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which Jg q = 0. Typically, one expects that 

Mnnc\ N=0 = A c = 2286.46 ± 0.14 MeV, (45) 
Mnnb\ N=0 = A b = 5620.2 ± 1.6 MeV. (46) 

These values are reproduced by formula (144]) with m c = 1.252 GeV and mj = 4.612 GeV. 
These masses, obtained by a comparison of the quark model to the experimental data, are 
clearly compatible with those obtained from the mass combination fill I) - both determina- 
tions actually differ by less than 5%. This is a first evidence of the compatibility between 
quark model and large N c expansion in the heavy baryon sector. The supplementary pa- 
rameters for qqQ baryons are gathered in Table IVIIII One can notice that we predict 
M nsc \ N=Q = 2433 MeV and M ns b\ N=0 = 5767 MeV with these parameters. These values are 
very close to the experimentally observed masses of H c and E b respectively. 



TABLE VIII: Supplementary parameters for qqQ baryons. 



Fixed parameters 


Fitted parameters 


fei = 0.930 


m c = 1.252 GeV 


a\ = 0.7ao 


m b = 4.612 GeV 



V. COMPARISON OF THE TWO APPROACHES 

First we recall that the heavy quark masses can be obtained by two different ways. On 
the one hand, the large V c inspired mass combination ffTTl) leads to m c = 1315 MeV and 
rrib = 4642 MeV. On the other hand, the quark model mass formula (l4*TT) is compatible with 
the experimental data provided m c = 1252 MeV and m b = 4612 MeV. Both approaches 
lead to quark masses that differ by less than 5%, as pointed out in Sec. IIVFI Thus the two 
approaches that are considered in this paper agree at least at the dominant order, where 
only itlq is present. 

The other parameter involved in the large N c mass formula is A. A comparison of the 
spin independent part of the mass formulas (J5j) and ( fldl) leads to the following identification 
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where m = ^3k x aj2. According to Eqs. © and (fT2bl) one has c = A ~ 0.324 GeV. The 
quark model parameters of Tables [VTT1 and [Villi give 0.333 GeV for the expression after the 
second equality sign in Eq. (j47|) . which means a very good agreement for the QCD scale 
A. In this quantity, 0.475 GeV comes from the dynamics of the confinement (4/^/3), while 
the Coulomb interaction (term containing ao and ai) contributes with —0.044 GeV and the 
self-energy (term proportional to /) with —0.097 GeV. The mass shift yielded by these two 
residual interactions is quite significant and their presence improves the value of A. 
Next the terms of order 1/rriQ lead to the identity 
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(48) 



2nn\) 6 V" J ' 3 
Note that to test this relation the value of rriQ is not needed, like for the identity (T47j). The 



large N c parameter A = 0.324 GeV gives for the left hand side of 



A 2 = 0.096 GeV 2 



and the quark model gives for the right hand side 0.091 GeV 2 , which is again a good 
agreement. In this quantity, the contribution of the dynamics of the confinement [fciaF(O)] 
is 0.105 GeV 2 , while the contributions of the Coulomb interaction and of the self-energy 
are 0.029 GeV 2 and —0.043 GeV 2 respectively. The relative magnitude of these two terms 
compared to the first one is here larger but they nearly cancel each other. 

The SU(3)-flavor breaking term is proportional to the factor eA x ~ m s — m in the 1/N C 
mass formula f)13p and similarly by the nonvanishing strange quark mass in the quark model 
where one takes m — 0. Using Eqs. ( fTBl . ( fT4l) . and ( |4T1) one obtains 



2m 2 
V3^i 
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12/ii 



ha f r- \ fa ( 3 (3 



(49) 



From phenomenology, Eq. (1141) implies that eA x = 0.206 GeV and the quark model estimate 
is 0.170 GeV, which compares satisfactorily with the value used in the combined 1/A^ C 
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and 1/rriQ expansion 12j. In the quark model, the contribution of the dynamics of the 
confinement (term proportional to 1/2) is 0.093 GeV, while the contributions of the Coulomb 
interaction and of the self-energy are —0.009 GeV and 0.085 GeV respectively. Thus the 
effect of the self-energy is as large as that of the confinement. 

Let us recall that, except for m c and nib, all the model parameters are determined from 
theoretical arguments combined with phenomenology, or are fitted on light baryon masses. 
The comparison of our results with the 1/N C expansion coefficients Cq, c and eA x are inde- 
pendent of the rriQ values. So we can say that this analysis is parameter free. 

So far, our formalism is spin- independent. An evaluation of the coefficients c-i, c' 2 , and 
C2 through a computation of the spin-dependent effects within a three-body quark model 
is then de facto out of the scope of the present approach. If included, the spin-dependent 
interactions between quarks i and j would appear as relativistic corrections to the Coulomb 
potential, proportional to l//ij/ij (fii is the dynamical mass of the quark i) [5]. At the 
dominant order, one expects that C2 oc [i{ 2 and c'^ oc fii 1 . The ratio 0^/02 should thus be 
of order = 356 MeV, which is roughly in agreement with Eq. ([9]) stating that d^jc^ ~ A. 
This gives an indication that the quark model and the 1/N C expansion method would remain 
compatible if the spin- dependent effects were included, as we already pointed out in the light 
baryon sector 5), Q] . 

Charm and bottom baryons have been studied with a Hamiltonian similar to ours in 



Ref. [31]. All parameters (a, as, m n , etc.) have values very close to ours but some differ- 
ences exist: A genuine junction Y is used for the confinement instead of the approximation 
(117p - (ll8j) . the auxiliary fields are introduced only at the level of the kinetic part, the Coulomb 
potential is not treated perturbatively and the color-magnetic interaction is taken into ac- 
count. The consequences of this procedure is that no analytical mass formula can be derived 
explicitly. But, the numerical results obtained in that paper are in good agreement with 
experiment, which reinforce our approach. Moreover it was also found that a unit of angular 
momentum between the heavy quark and the two light quarks is energetically favored with 
respect to a unit of angular momentum between the two light quarks. This correspond to 
our choice = 0. 
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VI. CONCLUSIONS 



Our previous studies establishing a connection between the quark model and the 1/N C 
expansion for light baryons have been successfully extended to baryons containing a heavy 
quark. Accordingly, the 1/N C expansion was supplemented by an 1/ tuq expansion due to the 
heavy quark. As in the light baryon sector, there is a clear correspondence between various 
terms appearing in our mass formula (|41|) and those of the mass formula in the combined 
1/N C and I/tuq expansion described in Sec. III. First, both methods lead to compatible 
values for the heavy quark masses. Second, the typical QCD scale involved in the 1/N C 
expansion is well reproduced by the quark model without any free parameter: All necessary 
parameters have been previously fitted on light baryons. Finally, the dominant term in 
SU(3)-flavor breaking expansion is satisfactorily reproduced. The spin-dependent terms, 
seen as relativistic effects, deserve a special study, to be considered in the future. 

This study, completing the two previous ones [5|, |6| , brings reliable QCD-based support in 
favor of the constituent quark model assumptions due to the compatibility of its mass formula 
and the mass formula derived from the model independent 1/N C expansion. Moreover, better 
insight into the coefficients q encoding the QCD dynamics in the mass operator is obtained: 
the dependence on the quark content and on the excitation number. 

We presently focused on ground state heavy baryons. For excited states, the quark 
model suggests that the band number N classifying the heavy baryon resonances should be 
associated to the quantum of excitation of the heavy quark-light diquark pair in a harmonic 
oscillator picture. We leave a detailed study of excited heavy baryons for future studies. 
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